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Magnetism of gadolinium: a first-principles perspective 
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By calculating the spectral density of states in the ferromagnetic (FM) ground state and in the 
high temperature paramagnetic (PM) phase we provide the first concise study of finite temperature 
effects on the electronic structure of the bulk and the surface of gadolinium metal. The variation 
of calculated spectral properties of the Fermi surface and the density of states in the bulk and 
at the surface are in good agreement with recent photoemission experiments performed in both 
ferromagnetic and paramagnetic phases. In the paramagnetic state we find vanishing spin splitting 
of the conduction band, but finite local spin moments both in bulk and at the surface. We clearly 
demonstrate that the formation of these local spin moments in the conduction band is due to the 
asymmetry of the density of states in the two spin channels, suggesting a complex, non-Stoner 
behavior. We, therefore, suggest that the vanishing or nearly vanishing spin splitting of spectral 
features cannot be used as an indicator for Stoner-like magnetism. 


Pure hep Gd metal is perhaps the most investigated 
strongly correlated metallic system, where the well lo¬ 
calized magnetic moment of the 4/-shell interacts with 
the itinerant electrons of the conduction band, jl], [3] The 
half-filled 4/-shell of Gd possesses a spin-moment of 7 hb, 
and it is energetically well separated from the conduction 
band of the spd-electrons. In the FM ground state the 
ordered local moments induce a spin splitting of about 
1 eV of the conduction band resulting in an additional 
spin magnetic moment of 0.6 f^B Q- 

Modern electronic structure theory describes ade¬ 
quately the ground state of Gd, if the strong correlation 
in the 4/-shell is treated in some, even simple, manner. 
In pioneering works it has been shown that the LSDA+U 
method provides a good description of the Gd ground 
state in bulk Hi and on the (0001) surface [6]. Later 
studies (j, 0 have confirmed this conclusion. Moreover, 
recent angle-resolved photoemission (ARPES) measure¬ 
ments of the Fermi surface in the FM phase of the bulk 
0 have been successfully reproduced using the LSDA+U 
methodology[l0]. Despite more sophisticated treatments 
of the correlated 4/-shell in Gd ]11-Il3j. the considerable 
splitting of the Gd conduction band could also be ob¬ 
tained in a remarkably simple spin-polarized open-core 
approach, [l^ 

Owing to the strong localization of 4/-states, the finite 
temperature magnetism in hep Gd seems to be well de¬ 
scribed by the Heisenberg model, and it may be regarded 
as a model Heisenberg system among metallic ferromag- 
nets. The magnetic critical temperature of hep Gd can 
indeed be estimated reasonably well in terms of a Heisen¬ 
berg model and calculating the exchange constants using 
either the open-core fl5U17l ]. the LSDA+U [10 or the 
self-interaction correction (SIC) approach [10 . 

However, despite the successful understanding of the 
ground state and finite temperature magnetism of Gd, 
the fundamental issue concerning the interaction between 


localized 4/-moments and itinerant electrons has been 
the subject of a heated debat e. fill . 1^-22| Earlier pho¬ 
toemission (PES) experiments [23j have shown that the 
spin splitting of the conduction band vanishes near the 
Curie temperature. This gives rise to a simple inter¬ 
pretation that the spd-bands experience a spin polar¬ 
ization due to the (average) magnetization of the lo¬ 
calized 4/-electrons, a hallmark of Stoner magnetism 
23|- 271. However, more precise analyses of later PES 
measurements [28l430( seemed to provide clear evidence 
that the spin splitting of majority and minority spin 
channels remains finite in the PM state, although very 
small, being on the verge of the PES resolution. From 
this observation it was concluded that the Stoner model 
cannot describe the magnetism of bulk Gd.f+O] 

Contradictory theoretical interpretations were also 
given based on ab initio determined model parameters 
for the description of the interaction between the local¬ 
ized 4/-moments and the conduction band. For example 
the authors of Ref. (31] derived a spin-mixing behavior of 
the conduction band in agreement with the experimental 
conclusions by Maiti et al. 21, 20, whereas later on the 


same authors proposed a theory predicting Stoner be¬ 
havior 0 in line with old measurements by Kim et al. 
[Til] , Fully ab initio attempts to resolve the above the¬ 
oretical contradiction used the disordered local moment 
approximation (DLM) to model finite temperature mag¬ 
netic disorder [221 ] or employed a study of non-collinear 
spin-configurations. [3] Although both of these studies 
were based on the open-core treatment of the 4/-shell, 
they indicated a spin-mixing behavior of the conduction 
band. The DLM calculation [H] predicted a small but 
finite splitting of the majority and minority spin chan¬ 
nels and a corresponding finite value of the local moment 
due to the spd-band in the paramagnetic state. The non- 
collinear calculation 0 clearly showed that the conduc¬ 
tion band moment is strongly coupled to the direction of 
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the local 4/-moment and cannot be treated as an inde¬ 
pendent magnetic degree of freedom. 

In recent years the study of ultrafast magnetic pro¬ 
cesses induced by femtosecond laser pulses has been fuel¬ 
ing a reinvigorated discussion concerning the nature of 
the interaction between 4/-states and the conduction 
band moments in pure Gd and its compounds. [32 — 


|34J A crucial point in the simulations of these processes 
lies in the assumption of treating 5 d- and 4 /-m oments 
as independent magnetic degrees of freedom. [351, [36J An 
elaborated discussion of this topic was recently given by 
Sandratskii.|f3| 

Recent ARPES 371 and positron annihilation spec¬ 
troscopy [11115 experiments were able to observe differ¬ 
ences of the Fermi surface topology between the FM and 
PM state. However, the current theoretical understand¬ 
ing of the PM electronic structure of the Gd conduction 
bands is far from being satisfactory. The Fermi surface 
of pure bulk Gd in the PM state has been calculated 
for several lattice constants using self-interaction correc¬ 
tion to the LSDA llj. The main features of the derived 
PM Fermi surface seem to be consistent with the later 
measured one [37|. However, so far there has been no 
consistent theoretical description of the Fermi surfaces 
in the FM and PM state and the evolution of spectro¬ 
scopic properties from bulk to surface within the same 
framework. 

An equivalently important issue is that the main exper¬ 
imental methods that provide us with precise information 
concerning the electronic structure of the Gd conduction 
band are surface sensitive. The presence of the surface 
leads to essential changes in the electronic structure of 
the terminating layers compared to the bulk. Moreover, 
the spin splitting of the surface states has been regarded 
as an important source of the information concerning the 
magnetic behavior of the conduction electrons. [19] It is 
therefore demanding to disentangle the signatures of the 
surface state from those of the bulk. 

In this Letter we address both issues, namely the 
changes of the electronic structure of hep Gd induced by 
finite temperature magnetic disorder and the variation of 
the Fermi surface near the Gd (0001) surface termination. 
We employ the well established LSDA+U methodology 
for the treatment of the half filled 4/-shell and calculate 
bulk and layer resolved Fermi surfaces, as well as densities 
of states in the FM and the DLM state. We show that 
the high temperature behavior of the conduction band is 
governed by spin-mixing both in the bulk and near the 
surface. The comparison of our results with PES 


21 , 2 ? 


and ARPES studies suggests a good agreement for 

both FM and PM phases with regard to bulk and surface 
properties. We also find that electronic states associated 
with a presence of the surface decay faster towards the 
bulk in the paramagnetic phase than in the ferromagnetic 
state. 

We performed self-consistent calculations by using 


the fully relativistic screened Korringa-Kohn-Rostoker 
(KKR) Green’s function method with the atomic sphere 
approximation, that allows the study of layered systems 
and surfaces. [41[ The strong correlation of the local¬ 
ized 4/-states was treated within the framework of the 
LSDA+U approach[3, fe] as implemented within the KKR 
method. [45 The paramagnetic phase is treated in terms 
of the relativistic disordered local moment (R-DLM) 
method. [43l. H] in which the coherent potential approx¬ 
imation (CPA) is employed to model disorder in local 
spin orientations. [45j The application of these meth¬ 
ods allows us to treat the correlations in the 4/-shell on 
equal footing in FM ground state and finite temperature 
PM phase. The spectral density of states (SDOS) in the 
FM phase is directly related to the KKR Green’s func¬ 
tion, while in the PM phase the SDOS is evaluated from 
the configurationally averaged Green’s function [4~0. ITgI|. 
For all the calculations we considered a hexagonal closed 
packed lattice with the experimental value of the c/a ra¬ 
tio of 1.5904, while the lattice constant was optimized to 
a = 3.450 A. The optimization was performed with the 
commonly used U = 6.7 eV and J = 0.7 eV values, [ij 
and the LSDA parametrization from Ref. [47| . 



Figure 1. a) Spin resolved spectral DOS of the ferromagnetic 
phase and b) total spectral DOS of the paramagnetic phase of 
bulk gadolinium at the Fermi energy as projected to the 2D 
Brillouin zone. Darker colors represent larger values of the 
spectral DOS. 

The spin resolved SDOS of the FM bulk and the total 
SDOS of the PM bulk is depicted in Fig. [Q projected 
onto the hexagonal two-dimensional (2D) Brillouin zone. 
As can be clearly seen in Fig. |T] a), the Fermi surface 
in the case of the FM phase shows a significant spin 
splitting in momentum space with concentric features 
around both the T and K points of the Brillouin zone, in 
agreement with ARPES experiments.[37j The minority 
(j.) components form a tight cylinder around the T point 
surrounded by a second cylinder consisting entirely of 
majority (f) spin-states. A similar, but much narrower 
feature can be observed around the K points with the 
order of the spin channels reversed, that is the majority 
states form the inner and the minority components the 
outer cylinder. In the paramagnetic phase the splitting 
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Figure 2. Spin resolved densities of states in Gd bulk and in 
the top layer of Gd(0001) surface in the a) ferromagnetic and 
b) paramagnetic phase. 


disappears, only a single cylindrical feature is present 
around the T and K points, nicely recovering the experi¬ 
mental observations.[37] These experimental and theoret¬ 
ical findings seem to lead to the simple conclusion that a 
vanishing spin splitting in the paramagnetic phase is an 
indication of Stoner behavior, where the spin splitting 
depends on the total magnetization of Gd. 

The densities of states for the bulk and surface in the 
FM and PM phases are depicted in Fig. [21 In the FM 
phase, similarly to the results of previous works 
the 4/-electrons move away from the Fermi level with an 
exchange splitting of ~ 11 eV. Note that the relatively 
large dispersion of about 0.8 eV of the majority 4/-states 
is due to spin-orbit splitting of these states. Clearly, in 
the PM phase very similar features of the density of 4/- 
states can be found. The conduction electrons, domi¬ 
nated by 5d-states, are characterized by a spin splitting 
on the order of 1 eV in the FM state. The magnitude 
of this splitting can be inferred from Fig. [2Ji), for exam¬ 
ple by comparing the onset of the majority d electrons 
at roughly — 3 eV below the Fermi energy to the onset of 
minority d-electrons at — 2eV. Another measure of the 
spin splitting would be to consider the shift of the min¬ 
imum in the majority-spin DOS at the Fermi energy to 
roughly 1.2 eV above the Fermi level in the minority-spin 
DOS. The spin splitting of bulk states in the PM phase 
is somewhat harder to quantify with the above measures. 
A comparison of the position of the peak near the Fermi 
energy for the two spin channels suggests a largely re¬ 
duced splitting of the order of 0.1 eV, as can be assessed 
from Fig. [2] b). At this point one might again conclude 
that an almost vanishing spin splitting in the PM phase 
indicates Stoner behavior. 

In the FM phase we calculated a total spin moment of 
7.77 mb, that is, the polarization of the spd -band amounts 


to 0.77 mb- In the PM phase the total spin magnetic mo¬ 
ment is reduced to 7.41 mb, which means that a consid¬ 
erable polarization of the conduction band still persists. 
These values are in agreement with recent disordered 
local moment calculations [ij], [22j and also with studies 
performed for non-collinear magnetic configurations 14|, 
where the 4/-electrons are treated as part of the core. 
The finite value of the local moment related to the con¬ 
duction band in the PM phase is indeed at the heart of 
the controversy, since there is (nearly) no splitting in the 
spin channels either in momentum space or in the energy 
resolved spectra. 

The density of states at the surface, shown in Fig. [21 
has two main characteristics. First, in accordance with 
previous results [6] and experimental observation, j48j the 
localized 4/-states experience a down-shift in energy. 
This feature is slightly more pronounced in the PM phase 
than in the FM phase. The second feature, again in line 
with experiments, [2lj, [28|, [30| is the appearance of new 
states absent in the bulk. The surface states in the FM 
phase near the Fermi energy are highly of majority-spin 
character. In the PM phase, we still detect a track of such 
surface states near the Fermi energy, but with a reduced 
spin polarization as observed in experiments. (2l[|2^| The 
spin magnetic moment per atom on the surface in the 
FM phase is found to be slightly reduced, compared to 
the bulk value, to 7.74 mb- On the other hand, in the PM 
phase the magnetization slightly increases to 7.48 mb on 
the surface. Therefore, similarly to the case of bulk, the 
vanishing splitting of the spectral density of the surface 
states in the PM phase can not be regarded as a herald 
of Stoner magnetism. It is the asymmetry of the spectral 
weight for the two spin channels that is responsible for 
the finite residual local moment of the conduction band 
in the PM phase, suggesting a spin-mixing non-Stoner 
behavior. 

As we stated before, most spectroscopic methods are 
surface sensitive. So far it is not clear how deep the 
effects of surface termination on the electronic spectra 
extend in to the bulk. To elucidate an answer for this 
question, in Fig. [3] we present the layer resolved spectral 
DOS at the Fermi energy for the first six layers below the 
(0001) surface and compare them with the bulk spectral 
function. The most pronounced effects of the surface can 
be seen on the first (surface) layer. In both phases a 
high intensity feature is present in the vicinity of the T 
point which is absent in the subsequent layers. In the FM 
phase an enhancement of the spectral density around the 
K points and along the lines connecting equivalent I\ 
points, forming a hexagram like structure, can also be 
attributed to the surface states. These highly spin polar¬ 
ized features gradually fade in deeper layers. In the PM 
phase no such distinctive spectral features are present 
near the surface. We merely observe a slight modulation 
in the SDOS along the T — K direction in the region out 
of the disk of zero weight characterizing the bulk. Our 
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Figure 3. Spectral densities of states evaluated at the Fermi energy in the ferromagnetic (top row) and paramagnetic (bottom 
row) phases for Gd bulk and the first six layers below a Gd(0001) surface. 


calculations show that the penetration depth of the sur¬ 
face states is different in the two phases. A bulk like 
spectral density is nearly recovered in the 3rd layer from 
the surface in the PM phase, while in the FM phase the 
SDOS even at the 5th-6th layers below the surface con¬ 
tain surface associated features. Thus one can deduce 
that the surface effects are more pronounced in the FM 
phase than in the high temperature PM state. 

In conclusion, we provided a coherent first-principles 
study of the electronic and magnetic properties of the 
bulk and the surface of Gd metal, both in the ferromag¬ 
netic ground state and in the high temperature param¬ 
agnetic phase. Our results reveal the root of the existing 
controversy concerning the nature of magnetism related 
to the conduction band of Gd. A finite, though reduced, 
value of the spin moment related to the conduction band 
in the paramagnetic phase both in the bulk and for the 
surface implies a non-Stoner magnetism for these elec¬ 
trons. We argue that the asymmetry of the density of 
states with respect to different spin channels can be re¬ 
garded as an indicator for this behavior rather than the 
spin splitting of the bulk conduction band or of the sur¬ 
face states, as it has been judged in many previous inves¬ 
tigations. The asymmetry of the conduction spin-bands 
is generated by stable local 4/-moments irrespective of 
their directions and it is the consequence of the forma¬ 
tion of a common band, i.e. of a spin-mixing state. 
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